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Abstract. This paper studies the combinatorial Yamabe flow on hyperbolic 
surfaces with boundary. It is proved by applying a variational principle that 
the length of boundary components is uniquely determined by the combina- 
torial conformal factor. The combinatorial Yamabe flow is a gradient flow of 
a concave function. The long time behavior of the flow and the geometric 
meaning is investigated. 



1. Introduction 

1.1. Piecewise flat metrics. In trying to develop the analogous piecewise linear 
conformal geometry, Luo studied the combinatorial Yamabe problem for piecewise 
flat metrics on triangulated surfaces jLuo04j . We summarize a part of this work in 
the following. Suppose S is a connected orientable closed surface with a triangula- 
tion T so that V,E,F are sets of all vertexes, edges and triangles in T. We identify 
a vertex of T with an index, an edge of T with a pair of indexes and a triangle of 
T with a triple of indexes. This means V = {1,2, ...n}, E = {ij \ i,j G V} and 
F = {ijk I k e V}, where n is the number of vertexes. 

A piecewise flat metric on (E, T) is identified with a vector indexed by the set of 
edges E. More precisely, it is an assignment to each edge a positive number such 
that the three numbers assigned to the three edges of a triangle satisfy the triangle 
inequality. Equipped with a piecewise flat metric, each triangle of T can be realized 
as a Euclidean triangle and E can be realized as a Euclidean polyhedral surface. 

Let's fix a piecewise flat metric on (E,T) as f' e IR.|^- The assignment to the 
edge ij is denoted by l^y A combinatorial conformal factor on (E,T) is a vector 
w = {wi,W2, ■■■,Wn) G R" which assigns each vertex i £ V a, number Wi. (In |Luo04j . 
the notation Ui is used, where Ui = e*"'.) From a combinatorial conformal factor, 

I E I 

we obtained a new vector I £ MI^q as follows: 
(1) Z,,=e-'+-^/?, 

for each edge ij S E. 

Let We be the space of combinatorial conformal factors w such that each vector I 
corresponding to a vector w € We is indeed a piecewise flat metric. In other words, 
the triangle inequality holds for each triangle under the assignment I. Obviously, 
We depends on the initial metric 

For a vector w S We, the vector I corresponding to w is a piecewise flat metric. 
Each triangle of T is realized as a Euclidean triangle under the metric / . At a vertex 
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i, the curvature of the metric I is defined as follows. Let be the inner angle of 
the triangle ijk € F between the edges ij and ik. Then 




sending a combinatorial conformal factor to the curvature. 
Theorem 1 (Luc). The map ipE is a local homeomorphism. 

The theorem is proved by applying a variational principle. An local convex 
energy function is constructed using the derivative cosine law. And V'e truns out 
to be the gradient of the energy function. 

Motivated by establishing a discrete Uniformization Theorem, Luo introduced 
the combinatorial Yamabe flow 



Corollary 2 (Luo). The combinatorial Yamabe flow (0j is the negative gradient 
flow of a locally convex function in terms of w. And J27=i decreasing in 

time t. 

1.2. Related work. Motivated by the application in computer graphics, Spring- 
born, Schroder, and Pinkall jSSPOSj considered this combinatorial conformal change 
of piecewise flat metrics ([T]). They found an explicit formula of the energy func- 
tion. Glickenstein [GliOSaj IGliOSb] studied the combinatorial Yamabe flow on 3- 
dimensional piecewise flat manifolds relating to the ball packing metric of Cooper 
and Rivin |CR96j . Recently Glickenstein |Gli09| set the theory of combinatorial 
Yamabe flow of piecewise flat metric in a broader context including the theory of 
circle packing on surfaces. This combinatorial conformal change of metrics has ap- 
peared in physic literature [RW84 and numerical analysis literature [Kcr96, PC98]. 
We were informed by Luo in 2009 that Springborn considered the combinatorial 
conformal change of hyperbolic metric on a triangulated closed surface. He intro- 
duced the combinatorial conformal change as 



1.3. Hyperbolic metrics. In this paper we study the combinatorial Yamabe flow 
on hyperbolic surfaces with geodesic boundary. Let S be a connected orientable 
compact surface with n boundary components. The set of boundary components 
is _B = {1, 2, ...n} where a boundary component is identified with an index. 

A colored hexagon is a hexagon with three non-pairwise adjacent edges labeled by 
red and the opposite edges labeled by black. Take a finite disjoint union of colored 
hexagons and identify all red edges in pairs by homeomorphisms. The quotient is 
a compact surface with non-empty boundary together with an ideal triangulation. 
The faces in the ideal triangulation are quotients of the hexagons. The quotients 
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of red edges are called the edges of the ideal triangulation while the quotients of 
black edges are called the boundary arcs. 

It is well-known that each connected orientable compact surface E of non-empty 
boundary and negative Euler characteristic admits an ideal triangulation. 

Let T be an ideal triangulation of S. Since an edge of T connects two boundary 
components of S, an edge of T is indentified with a pair of indexes. The set of edges 
of T is -E = {ij I i, J G i?}. Since a face in T is determined by its boundary arcs in 
three boundary components, a face of T is indentified with a triple of indexes. The 
set of faces oiT \s F = {ijk \ k G B}. 

A hyperbolic metric on (S,T) is identified with a vector indexed by the set of 
edges E. More precisely, it is an assignment to each edge a positive number. It is 
well-known that |Bus92', for any three positive numbers, there exists a hyperbolic 
right-angled hexagon the length of whose three non-pairwise adjacent edges are the 
three numbers. Furthermore, the hexagon is unique up to isometry. Therefore, for 

I -El 

vector in M!^q , each face of F can be realized as a unique hyperbolic right-angled 
hexagon and the surface E can be realized as a hyperbolic surface with geodesic 
boundary. 

Let's fix a hyperbolic metric on (E,r) as P G I^m)- The assignment to the 
edge ij is denoted by l^y A combinatorial conformal factor on (E,T) is a vector 
w = (wi, W2, Wn) G which assigns each boundary component i G -B a number 
Wi. From a combinatorial conformal factor, we obtained a new assignment I G M'^ 
as follows 

(4) cosh ^ = e"-+'"^ cosh ^ 

for any edge ij G E. This definition (|4]) is an analogue of Springborn's definition 
of combinatorial conformal change of hyperbolic metrics on triangulated closed 
surfaces (jS]). 

Denote by W the set of combinatorial conformal factors such that the corre- 
sponding assignment is positive, i.e., I G R^. Therefore, for w G W, we obtained 
a new hyperbolic metric on (E,r). Each face of T is a hyperbolic right-angled 
hexagon. Let 9jf, be the length of the boundary arc in the boundary component i 
of the face ijk G F. Then the length of the boundary component i is 

This produces a map 

sending a combinatorial conformal factor to the length of boundary components. 

Theorem 3. The map ip is a homeomorphism. 

This is a result of global rigidity while Theorem 1 is a result of local rigidity. 
We also consider the combinatorial Yamabe fiow in this situation 
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Corollary 4. The combinatorial Yamabe flow 0) is the gradient flow of a concave 
function in terms of w. And ^ii^) decreasing in time t. 

We investigate the long time behavior of the flow. 

Theorem 5. The combinatorial Yamabe flow ^ has a solution for t £ [0,oo). 
Along the flow any initial hyperbolic surface with geodesic boundary converges 
to a complete hyperbolic surface with cusps. 

f.4. Variational principle. The approach of variational principle of studying 
polyhedral surfaces was introduced by Colin de Verdiere |CdV91) in his proof 
of Andreev-Thurston's circle packing theorem. Since then, many works about 
variational principles on polyhedral surfaces have appeared. For example, see 
|Bra92l iRwMl iLelOl ICLfM IBHImI IHIoOfil IGuo07l [Sp?58l [GLOQI IGuoOQ) and others. 

1.5. Organization of the paper. Theorem[3]is proved in section 2. Corollary |4] 
and Theorem [5] are proved in section 3. 



2. HOMEOMORPHISM 

2.1. Space of combinatorial conformal factors. Let l^ S k|^' be a fixed hy- 
perbolic metric on (E,r). We investigate the space of combinatorial conformal 
factor such that the corresponding new assignment is a hyperbolic metric. For a 
face ijk e F, denote by W*^*^ the space of vectors {wi,Wj,Wk) such that ljk,lki and 
lij are positive. 

Lemma 6. W'-''^ is a convex polytope. 

Proof. By definition Q cosh = e^''^^^ cosh The only requirement is kj > 0. 
Hence 

Wi + W j > — In cosh — . 

2 

Similar inequalities hold for wj+Wk and Wk + Wi. Therefore W'-''^ is the intersection 
of three half space. □ 

Corollary 7. The space W is a convex polytope. 

Proof yy = njjfeGFW'^''- □ 

2.2. Energy function. Let's focus on one face ijk £ F. When {wi,Wj, Wk) £ W*-''"', 
from (|4]), we obtain a hyperbolic right-angled hexagon whose non-pairwise adjacent 
edges have length Ijk, hi and kj. Let 0]k^^ki ^^"^ ^ij length of the hyperbolic 
arcs opposite to the edges jfc, ki and ij of this hexagon. Therefore dj^,6l,- and 6fj 
are functions of Wi , Wj , Wk . 

Lemma 8. The Jacobian matrix of functions d^jk,(^ki'^ij *^ terms of Wi,Wj ,Wk is 
symmetric. 
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Proof. The cosine law for hyperbolic right-angled hexagon induces the derivative 
cosine law: 

dOlk \ _i / sinh/,fc 

del = 1 sinh/fe, 

4J ) sinh^?^.sinhZfc,sinh/,fc \^ q ^ ^.^^^^ 

-1 cosh^fe cosh^^^ \ / dl,k 
cosh0*j —1 cosh 6*^ J. I dlki 
cosh6'(.^ cosh6'ifc -1 /V dUj 



By differentiating the two sides of the equation cosh ^ = e"''+"'j cosh-j-, 
we obtain 

2 sinhZi," 

"'ii = n Tidwi + dwj). 

cosh hj -1 

Similar formulas hold for dljk and dlki- Then we have 



dO]k 
d9% 




cosh 0, 



cosh - 



For simplicity of the notations, the above formula is written as 
dQ\ 





sinh 0* sinh l^i sinh / 



where M is a product of four matrixes. To prove the lemma, it is enough to show 
that the matrix M is symmetric. 

Represent cosh 0*^ , cosh 0^ - , cosh Q\- as functions of cosh 1^^ , cosh l^i , cosh lij using 
the cosine law. For simplicity of the notations, let a :~ coshl jk,b :— coshlki,c :— 
cosh lij , Then we have 

a+b-c+1 a+c-6+1 



M 



ah b 
+ — 



- 1 
1 



1 

-c-b 



c 

c + ab 

a- 1 

b + c- 



1 

a 



be 



\ 



1 



c - 1 
a+1 



a- 1 
ac a 



be 



□ 



Lemma 9. The Jacobian matrix of funetions &]kT^kv^tj terms of Wi,Wj ,Wk is 
negative definite. 
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Proof. Wc only need to show that the matrix M is positive definite. Let Mrs be 
the entry of M at r— th row and s— th column. First, Mu > 0. Second, 

c + ab 

Mil - Mi2 = + a + 1 > 0, 

0—1 

M22 - M21 = ^-^^ + 6 + 1 > 0. 

a — 1 

Then M11M22 - M12M21 > 0. Third, 



\2 



det M = sinh Ijk sinh Iki sinh • (sinh 0^ j. sinh 6*;^,- sinh 

smhljk sinhZfc,j sinh/ij 



cosh/jfe — 1 cosh/fci — 1 coshZy — 1 



2 > 0. 

□ 



Corollary 10. The differential 1-form Ojj^dwi + 61,-dwj + 9^jdwk is closed on W^^'^ 
For any c £ W*-''^', the integral 

£{wi, Wj,Wk) = / {6]i,dwi + Olidwj + O'ljdwk) 



is a strictly concave function on W*-''^ satisfying 

^=0' —^0^ —^0^ 

dw, '''' dwj dwu 

Proof. The differential 1-form is closed due to Lemma |S1 Since W'J'^ is connected 
and simply connected due to Lemma [6l then the function £{wi,Wj,Wk) is well 
defined, i.e., independent the path of integration. By LemmalHl the Hessian matrix 
of £{wi, Wj, Wk) is negative definite. □ 

2.3. Homeomorphism. In this subsection we prove Theorem[3lThe following two 
lemmas are needed. 

The first one is well-known in analysis. 

Lemma 11. Suppose X is an open convex set in M.^ and f : X M. a smooth 
function. If the Hessian matrix of f is positive definite for all x G X , then the 
gradient V/ : X — )■ is a smooth embedding. 

Lemma 12. For a family of combinatorial conformal factor w'^™' G W, if 
limm_^oo 'f^fe™'' = 00 for some k, then Ivoim^oo = and the convergence is 

independent of the values of Ivoim^oo w'f^^ for r ^ k. 



Proof. By definition (U), 



coshZj.™^ =e2"'r'+2"'^'ci-l, 



coshzi7Ue2-r'+2'"r'c2-i, 
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where ci — 2 cosh , C2 = 2 cosh -j- , C3 = 2 cosh -j- . Then 



k \{m) 
ij) 

cosh/}'"^ + coshZy/''' coshZt 



hm cosh(6'f ) 



1 . tj jk ki 

hm — ^ — 

sinh/]"^sinh/J^7^ 

..,,''"'4-91.,'™' , / 91.,'™' 4-9i.,'™' 9.,.(™' J^9,,.(™) 

hm 



;2-r'+2"'r'c2 - 2ci • e-^'+^^Ve'"^"' c2 - 2c2 



hm 



+ 1. 



(ci - e-2-<""-2-r')(c2 - e-2»<'"'-2»<-') 
+ hm — ^^^^^^^^^^^^^^=^^^^^^^^^^^^^^= 

/(c? - 2cie-2-y"'-^-^')(ci - 2c2e-2-r"'-2-i 



Hence hmm_+oo(^^ij)''™'' = independent the values of Yunm^a^Wr for r ^ k. 
Thus hm™^„o Si"'^ = lim™^oo EyfceF(^y )*™^ =0- ^ 

Proof of Theorem\^ Let /° G m|^' be a fixed hyperbohc metric on (E,r). For any 
combinatorial conformal factor w £ W, we obtain a new hyperbolic metric I e . 
By Corollary [ini for each face ijk e i^, there is a function £{wi,Wj ,Wk)- Define a 
function £ : W ^ R by 

f (Wl, U;2, Wn) = ^ £{wi,Wj,Wk) 
ijk£F 

where the sum is over all faces in F. By Corollary [101 £ is strictly concave on W 
and 

(6) E V e), = B.. 

That means the gradient of £ is exactly the map ip sending a combinatorial con- 
formal factor w to the corresponding length of boundary components. Thus ?/> is a 
smooth embedding due to Lemma [TTl 

To show that is a homeomorphism, we will prove that tpCW) is both open and 
closed in R"q. 

Since i/' is a smooth embedding, ipCW) is open in R"o- 

To show that ip{W) is closed in R"q, take a sequence of combinatorial conformal 
factor in W such that lim,„^oo(^J"\ -, si"^) e R^q. To prove the 

closeness, it is sufficient to show that there is a subsequence of w^™) whose limit is 
in W. 

Suppose otherwise, there is a subsequence, still denoted by tu^™', so that its 
limit is on the boundary of W. The first possibility is that there is some k such 
that limm_^oo Wi.™^ = 00. By Lemma ll2[ limm_j.oo S^™'* = 0. This contradicts the 
assumption that lim,„^oo(sJ'"', i?2™', ...,-Br™'') G ^>o- 
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The second possibility is that hmm^oo(w-"*'' + Wj'"^) — — In cosh ^ for some 
edge ij. That means hmm_j.oo 4]"' = 0- ■f'oi' the face ijk G F, we have 

. . s cosh /^.r^ + cosh /^r' cosh ^r"* 

™^oo m^oo sinh^r sinhZ.^J^^ 



> Um 



cosh ll^^ cosh /j^™^ 
sinh/^^"^sinhZ,^"^ 
cosh 

> hm = oo. 

sinh QT 

'■J 

Therefore \m\m-^oo{6]k)^™'^ — °° ^'^'^ hm„i^oo = oo. This contradicts the 
assumption that hm„^oo(-B^™\ 5^"^ si™^) e M'-^g. □ 

3. Flow 

Proof of Corollary [7} The combinatorial Yamabe flow ^ is the gradient flow of 
the concave function £{witW2, ...^Wn) due to the equation (j6|). 
Since 



^ dwi dt dwn dt dwk dt 

ijkeF ■' 

— ' dWi dwj dwk 



we have 



i=l 

^ — ' dWi dwj ■' dwk 



ijk£F 



.del dO]^^ dO' del, 

dWj dWk dWk dWi dWi dWj 

■ u^r, d[W^,Wj,Wk) 

By Lemma m the Jacobian matrix j*" 'j is negative definite. 

Hence ^4iT.7=i Bf{t) < 0. Therefore Bfit) is decreasing in t. □ 
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Proof of Theorem\^ Since B^{t) > 0, then Wi{t) > u;,(0) = 0. 
For any L < oo, we claim that liuit^L Wi{t) < oo. 

Otherwise, if liisit^L Wi{t) = oo for some L < oo, then by Lemma II 2[ we see 
liuit^L Bi{t) = 0. Therefore, for any e > 0, there exists some S > such that 
when t £ (L ~ e,L), the inequahties < Bi{t) < e holds. Hence, by the flow ([5]), 

< < e holds. Thus Wi{0) < Wi{t) < et < eL. This contradicts to the 
assumption that limt_j.i Wi{t) — co. 

Hence the solution of the flow ([5]) exists for all time t G [0, oo). 

To obtain the geometric picture, we claim that limt-^oo Bi{t) — for each 1 < 

1 < n. There are two cases to consider. 

First, if limt^oo Wi{t) = oo, by Lemma fT2l limt^oo Bi{t) = holds. 

Second, if limt_j.oo Wi{t) < oo, we claim limt_i.oc Bi{t) = still holds. Otherwise, 
Mmt^oo Bi{t) — a > 0. Then, for any e G (0,a), there exists some P > such 
that when t > P, the inequality Bi{t) > a ~ e holds. Hence, by the flow (O, 
> a — e holds. Therefore Wi{t) > {a — e)t. This contradicts to the assumption 
limt_^oo < oo. 

Once limj^oo Bi{t) = for each 1 < i < n, we have limt^oo (^]ki^) — for s-ny 
face ijk e F. Hence, for any edge ij, 



lim lij{t) — lim 



cosh0^,(t)+ cosh ej-.W cosh 



t-^oo--^^' t^oo sinh6'^,(t)sinh6l^^(t) 



J*. 

cosh g;, ft) cosh 0^,,(t) 

sinh0;,(^)sinh^^i,(^) 



> lim 



cosh ft) 
> lim — = oo. 

t^oo smh6'*^(i) 

Therefor each hyperbolic right-angled hexagon converges to a hyperbolic ideal tri- 
angle. □ 
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